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Abstract. The FENE dumbbell model consists of the incompressible Navier- 
Stokes equation and the Fokker-Planck equation for the polymer distribution. In 
such a model, the polymer elongation cannot exceed a limit Vb, yielding all interest- 
ing features near the boundary. In this paper we establish the local wcll-posedncss 
for the FENE dumbbell model under a class of Dirichlet-type boundary conditions 
dictated by the parameter b. As a result, for each b > we identify a sharp 
boundary requirement for the underlying density distribution, while the sharpness 
follows from the existence result for each specification of the boundary behavior. 
It is shown that the probability density governed by the Fokker-Planck equation 
approaches zero near boundary, necessarily faster than the distance function d for 
b > 2, faster than d\lnd\ for b = 2, and as fast as d h / 2 for < b < 2. Moreover, 
the sharp boundary requirement for b > 2 is also sufficient for the distribution to 
remain a probability density. 



1. Introduction 

Let N > 2 be an integer. We consider a dimer - an idealized polymer chain - as an 
elastic dumbbell consisting of two beads joined by a spring that can be modeled by 
an elongation vector m G M N (see e.g [6]), with \I/ being the elastic spring potential 
defined by 

Hb ( \m\ 2 \ 
(1.1) ^(m) = — — log ( 1 - j , meB. 

Here B : = 5(0, Vb) is a ball in ~R N with radius Vb denoting the maximum dumbbell 
extension. In the limiting case, this reduces to the Hookean model with \l/(m) = 
H\m\ 2 /2. A general bead-spring chain model may contain more than two beads 
coupled with elastic springs to represent a polymer chain. 

Polymers as such when put into an incompressible, viscous, isothermal Newtonian 
solvent are modeled by a system coupling the incompressible Navier-Stokes equation 
for the macroscopic velocity field v(t,x) with the Fokker-Planck equation for the 
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probability distribution function f(t,x,m) : 



(1.2) d t v + (v ■ V)v + Vp = 

(1.3) V-v = 



V ■ r + v k Av, 
0, 

-V m • (V m *(m)/) + 



(1.4) V)/ + V m -(VW) 



2fc fl r a 

C 



A m /, 



where x e K is the macroscopic Eulerian coordinate and m G B C IR is the 
microscopic molecular configuration variable. The model describes diluted solutions 
of polymeric liquids with noninteracting polymer chains (dimers). Note that the 
Fokker-Planck equation can be conveniently augmented to incorporate other effects 
such as inertial forces (see [T4]). 

In Navier-Stokes equation (II. 2p . p is hydrostatic pressure, is the kinematic vis- 
cosity coefficient, and r is a tensor representing the polymer contribution to stress, 



where X p is the polymer density constant. In the Fokker-Planck equation (11.41) . ( is 
the friction coefficient of the dumbbell beads, T a is the absolute temperature, and ks 
is the Boltzmann constant. We refer to [6j [131 H2] for a comprehensive survey of the 
physical background, and [42] for the computational aspect. 

Since B is a bounded domain, one has to add an appropriate boundary condition 
for / on the boundary dB. However, the singularity of the Fokker-Planck equation 
near dB makes the boundary issue rather subtle, and presents various challenges. 
To address the boundary issue, several transformations relating to the equilibrium 
solution have been introduced in literature ( see, e.g. [HJ [2U EH ESj) • It was shown 
in [33] that b = 2 is a threshold in the sense that for b > 2 any preassigned boundary 
value of the ratio of the distribution and the equilibrium will become redundant, and 
for b < 2 that value has to be a priori given. 

Our main quest in this paper is that what is the least boundary requirement for 
/ so that both existence and uniqueness of solutions to the FENE model can be 
established, also the solution remains a probability density. 

We addressed this issue in |36j for the microscopic FENE model alone and when 
b > 2. In this article we consider the well-posedness of the coupled system (11.21) - 
(jl.4p . A general discussion of this problem and background references are given in 
the introduction to [36J. Here we have two objectives: 

(1) to identify sharp boundary conditions on dB for all b > 0. 

(2) to prove well-posedness for the coupled FENE dumbbell model under the 
identified boundary condition. 

The setting for our problem is the coupled system subject to the initial data 





q(t,x,m)\ dB . 



v (x) 
fo(x,m), 
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Here v depends on b through the distance function, and q is a given function measuring 
the relative ratio of f /u near boundary. Our goal is to investigate solvability of the 
above system with the Cauchy-Dirichlet data. Note that our boundary condition is 
more or less a boundary behavior requirement for /, instead of the Dirichlet data in 
the traditional sense. 

Instead of using the distance function d = Vb — \m\ we shall use a regularized 
distance function p = b — \m\ 2 when describing the solution behavior near boundary. 
Our main observation is the form of v 



With some regularity requirement on q as well as on initial data we prove local well- 
posedness for the Cauchy-Direchlet problem in a weighted Sobolev space for each 
given q. Our results indicate that simply putting / = on boundary does not 
guarantee uniqueness of the solution. 

For the Dirichlet-type boundary condition (\1.7\i considered in this paper, our strat- 
egy is to study the transformed problem via 



with v defined in (II. 8ft so as to extract useful info for /. Inspired from |3H], for the 
coupled FENE system we use weak norm in m and strong norm in x, this enables us 
to prove wellposedness for all cases of b > and any given smooth q. 

For the case b > 2 of physical interest, we prove that / remains a density distribu- 
tion if and only if q\gs = 0. We thus identify a sharp boundary requirement for each 
b > for the underlying density distribution, while the sharpness is a consequence of 
the existence result for each q ^ 0. In particular, our result asserts that near bound- 
ary the probability density governed by the Fokker-Planck equation approaches zero, 
necessarily faster than the distance function d for b > 2, faster than d\lnd\ for 6 = 2, 
and as fast as d h l 2 for < b < 2. But within our current framework we have not been 
able to identify a non-trivial q for < b < 2 such that the corresponding solution is 
a density distribution. 

We remark that the sharp boundary condition presented in this work provides a 
threshold on the boundary requirement: subject to this condition or any stronger ones 
incorporated through a weighted function space [17] or just zero flux |3E], the Fokker- 
Planck dynamics will select the physically relevant solution, which is a probability 
density, any weaker boundary requirement can lead to many solutions, each depending 
on the ratio of f/u near boundary. 

This article is organized as follows. In Section 2, we state our main results and 
mains ideas of the proofs. In Section 3, we study the Fokker-Planck operator and 
well-posedness of the initial boundary value problem for the Fokker-Planck equation 
alone. This improves upon our previous work in [36]. The main result is summarized 
in Theorem [121 The Fokker-Planck problem involving spatial variable x is investi- 
gated in Section 4. Well-posedness of the coupled system is proved in Section 5. In 
Section 6, we sketch the proof of well-posedness for the coupled system with b > 6 in 
a different function space than what we used in Section 5. Some concluding remarks 



(1.8) 




< b < 2, 
6 = 2, 
b > 2. 



w = q 

v 
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are drawn in Section 7. 

We conclude this section by some bibliographical remarks. 

Existence results for the FENE model are usually limited to small-time existence 
and uniqueness of strong solutions. We refer to [H] for the local existence on some 
related coupled systems, [22] for the FENE model (in the setting where the Fokker- 
Planck equation is formulated by a stochastic differential equation) with b > 6, [17] 
for a polynomial force. More related to this paper are the work by Zhang and Zhang 
pET] for the FENE model when b > 76, and Masmoudi [38J for b > 0. Global existence 
results are usually limited to solutions near equilibrium, see [281 E2], or to some 
2D simplified models [TUJ d21 EZl E] • For results concerning the existence of weak 
solutions to the coupled FENE system we refer to [21 El SI El EH ESI SSI SH] ■ 

Boundary behavior of the polymer distribution governed by the FENE model is 
also essential in several other aspects, including the study of large time behavior, see 
[D HH [231 SS]; and development of numerical methods, see, e.g., [8J El M, [Ml EH M\ ■ 
We also refer to [21] for references on numerical aspects of polymeric fluid models. 

There are also some interesting works on non-Newtonian fluid models derived 
through a closure of the linear Fokker-Planck equation (see, e.g., [13 [IS])- We can 
refer to the pioneering work [181 [TUJ , and more recently to [lTJ [221 EDI ED I3"2| . 

However, none of these works is concerned with the sharpness of boundary condi- 
tions in terms of the elongation parameter. 



2. Main results 

After a suitable scaling and choice of parameters we arrive at the following Cauchy- 
Dirichlet probelm for the coupled system 

(2.1a) d t v + (v ■ V)v + Vp = 

(2.1b) V-v = 

(2.1c)d t f + (vV)f + V m -(Vvmf) = 

(2.1d) r = 

(2.1e) v{0,x) = 

(2-lf) /(0,x,m) = 

(2.1g) /(t,x,m)z/ _1 | ai? = 

To present our main results we first fix notations to be used throughout this article. 
We fix an exponent s, which is an integer in the range s > N/2 + 1. We use C to 
denote various constants depending on s, b and on some other quantities which we 
will indicate in the sequel. A 6-dependent weight function is defined as 

r p fe / 2 , < b < 2, 
(2.2) fjL=\ pln 2 f, 6 = 2, 

[ p 2 ~ h l\ b > 2. 



V-t + Av, t>Q, 
o, 

1 fbm \ 1 

-V m • (-/] + -A m /, meB, 

f bm , . 
/ m <8 — jam, 

J P 

fo(x,m), 
q(t,x,m)\dB- 



For b > 6, we also use 

(2.3) /i = p\ -1 < 9 < 1, b> 6. 

Other notations are listed as below as well. 

• L 2 ^ = |(/> : J (f) 2 ndm < oo j- 

. i^ = {0:0,\ i 0eLj,7 = l---iV.} 

o 

• if* denotes the completion of with if* norm. 

• if* is a dual space of H 

• if* is the usual Sobolev space with respect to x 



\v 



2 



\w 



0,s 



V / \d a v\ 2 dx, 

a\<s J ^ N 

/ / \d a w\ 2 fidmdx, 
~? Jrk Jb 



His = \ W \ls+ \VmU>\o,,> 

\\ w \\i,i,s = su P(klL + \ d Ml, 

t 

\\q\\ = Ikll^ + ll^gllifi. 

. ffJLj = {0 : |0| o , s < oo}, H° x Hl = {0 : |0| M < oo}. 

• L 2 H = L 2 ((0,T);if), C t H = C([0,T];H) for < t < T. 

• U = {0 : \ \0\\ L 2 H i + 110*11^^1), < oo} 

H = {0(t, •) 6 Hi : 1 10| | Lln + ll&ll^. < oo}. 

X, = [C t H s x n L 2 H^] x [C^LJ n L\H' a H%. 

• For a generic constant independent of T and a £ L 2 we denote 

(2.4) F(a) = c(r + J \a{t)\ 2 dt 

Due to such a constant, any two instances of F should be presumed to be with 
different constants. 

We now state our main theorem as follows: 

Theorem 1. Let b > and s be an integer such that s > N/2 + 1. Suppose that 
vo G H s x , /o^ 1 £ H^L 2 ^, and q e C}Hl +1 H^. Then, for some T > there exists a 
unique solution (v, f) to the coupled problem (12. ip such that 

o 

It is known from ([25]) that if* = if* for b > 6 with /x defined in (12.21) . Thus, 
the boundary condition Q2.1g ) is nothing but the zero dirichlet boundary condition 



under the assumption on q in Theorem [U For non-trivial q when b > 6, we show 
the well-posedness in a different weighted Sobolev space. The result summarized as 
below. 
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Theorem 2. Let b > 6 and s be an integer such that s > N/2 + 1. Suppose that 
v G H s x , f v~ l e H^L 2 ^, and q G C^H^H^ with fi defined in ( EJ[) . T/ien, /or 
some T > there exists a unique solution (v, f) to the coupled problem (12 .ip such 
that 

Theorem [1] and [2] tell us that for each given g, which denotes the rate of / approach- 
ing to zero relative to v near dB, there exists a unique solution (v,f). Also, they 
indicate that any weaker boundary requirement may lead to more than one solutions 
to (12.11) . For instance, the boundary condition 

iV~V £ Isb = 0, e > 

gives infinitely many solutions to (12. ip . Precisely we state the following non-uniqueness 
result. 

Theorem 3. Let v be a smooth function of p such that 
(2.5) lim- = 0. 



Then, the coupled problem (12. ip with ( |2.1g ) replaced by 



(2.6) f{t,x,m)i>- l \ d B = 

has infinitely many solutions in X M and X w for < b < 6 and b > 6 respectively. 

The natural question is for what q the obtained distribution / is a probability 
distribution. The answer when b > 2 is given in the following theorem. 

Theorem 4. Suppose that b > 2 and > 0. Under the assumption of Theorem 
U\ or\^ the unique solution f to the Cauchy-Direchlet problem (I2.ip is a probability 
distribution if and only if q\dB = 0. That is, f > if fo > 0, and for any t > 0, 

x e R N , 

(2.7) / f(t,x,m)dm= / f (x,m)dm. 

Jb Jb 

Theorem Q] is proven by a fixed point argument, which is now outlined. Given 

(u,g), we first solve the Navier-Stokes equation (NSE): 

(2.8a) d t v + (u ■ V)v + Vp = V • r + Av, 

(2.8b) V-v = 0, 

(2.8c) v(0,x) = v {x), 

f bm 

(2.8d) r = / m <g> — gdm. 



P 

With the obtained v we solve the Fokker-Planck equation (FPE): 

1 {bm \ 1 
(2.9a) 3t/+(T>-V)/ + V m -(Vwi/) = -V m ■ f — M + -A m /, 

(2.9b) /(0,x,m) = f (x,m), 

(2.9c) f(t,x,m)^ 1 \ dB = q(t,x,m)\ dB . 

The above two systems define a mapping (w,g) — > (v,f), the existence of problem 
(12. ip is equivalent to existence of a fixed point of this mapping. 



The main difficulty lies in monitoring the boundary behavior of /. Our strategy is 
to apply the transformation 

(2.10) f = u{w + q), 

to (12.91) to obtain a w-problem 

(2.11a) n(dt + v • V)w + L[w] = fih, 

(2.11b) w(0,x,m) = Wo(x,m), 

(2.11c) w(t,x,m)\gB = 0. 

Here the operator L is induced from the Fokker-Planck operator, v and /j, are weights 
depending on the distance functions defined in (jl.8p and (12. 2p respectively. The source 
term is obtained from q 

(2.12) h = -d t q-(v-V)q- f i,- 1 L[q], 
and the initial data is given by 

(2.13) Wo(x, m) := fo(x, m)v~ x — q(0, x, m). 
For given (u, w) with g = v(w + q), we arrive at a map T . 

J= : M -> M 
(u, w) i — y (v, w) 

Here M is a subset of 

C t H' x x [C t H°Ll n L?2£fTj] 

such that 

M = \(v,w): sup \v\ 2 s <Ax, sup \w\l s + l- [ \V m w\l s dt < A 2 \ . 

y 0<t<T 0<t<T * Jo J 

The strategy for the fixed point proof, which we implement in sections to follow, is to 
first prove that J 7 is well defined for some T, A\ and A 2 , then show that J 7 is actually 
a contraction map in a weak norm. Moreover, we will show that 

(2.14) J-(M) C X„. 
This proves Theorem [1] for 

q e ClH^Hl c [C t H°Ll n LtH*H*] . 
Theorem [2] is proved in the same manner. A sketch of proof is presented in Section 

6. 

In order to prove Theorem[3], we pick q(t, x, •) G C°°(B) f)C(B) and q\gB ^ such 
that 

C}H^Hl 0<b<6, 



n c. ) 1 x A" 

^ I C}H^H^ b > 6. 

Note that existence of such a g follows from the density of the weighted Sobolev 
space (see [25] for details). Then for each q we have a unique solution (y, f) to the 
coupled problem (12. lft from Theorem [1] and Theorem |2j Now, we check the boundary 
condition (12 .61) . 

\dB — JV —\QB — q-\dB, 
V V 
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which vanishes since q\g B is bounded and condition (12. 5 j) holds. This proves Theorem 

El 

Theorem H] follows from Proposition [02 and HH] via a flow map to be described in 
Section HJ The case for b > 6 can be proved by a simple modification, which is also 
sketched in Section 6. 

3. The Fokker-Planck operator 

We start with (12. 9p when x is not involved. In such a case it reduces to the following 
problem: 

(3.1a) d t f + C[f] = 0, meB,t>0, 

(3.1b) f(0,m) = / (m), 

(3.1c) f(t,m)i>~ l \ aB = q(t,m)\ dB . 

Here 

1 fbm \ 1 
(3-2) C[f] := V ■ (Kmf) - -V • (—/J - -A/, 

k = is a square integrable matrix function such that Tt(k) = 0. We omit m from 
V m in (I3.2p for notational convenience. 
The goal of this section is two folds; 

(1) to provide tools for subsequent sections. 

(2) to elaborate on this model alone as an extension of our previous work [36] . 
3.1. Transformed operator. The transformation (I2.10p leads to 

(3.3a) d t wfi + L[w] = /ih, m<EB,t>0, 

(3.3b) w(0,m) = wo, 

(3.3c) w(t,m)\ aB = 0, 

with the transformed operator L determined by 

(3.4) L[w] = iw^fijw]. 

The source term h = —d t q — [i~ l L[q\ and initial data for w is Wq = /o^ -1 — ?(0, m). 
From a direct calculation with the choice of \x in (I2.2p . and v in (11. 8p . (13. 4p can be 

expressed as 

(3.5) L[w] = --V • (Vwfi) + V • (Kmwfi) - Kw, 
where 

{0, 0<6<2, 
(JV + 2«m-m)ln-, 6 = 2, 

(N + 2Km-m)(b/2- l)p l ~ b/2 , b>2. 
Associated with the operator L, we define its time- dependent bilinear form 

(3.7) B[w,<f>;t] := J ^-Vw ■ — wpnm ■ V(f) — Kwcj^j dm 

o 

for 0, w G -H"^ and fixed t > 0. 
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We now describe the weak solution which we are looking for. 

o 

Definition 5. A function w EH is a weak solution of w-problem A3. 3\) . provided 

o 

(1) For each 4> G H 1 and almost every < t < T , 
(3.8) (d t w,<t>)o+B[w,<t>;t] = (h,<l>).. 

j-L (J. 

(2) w(0,m) = Wo(m) in L 2 sense, i.e. 

\w(0, m) — wo(m)\ 2 fxdm = 0. 



B 



Remark 6. In A3. 8\) . (ip,<f))o i is a dual pair for ip G (H X A* and cf) G H 1 , and can 

Hp 

be regarded as L 2 inner product. Indeed, from the Riesz representaiton theorem, for 

o o 

each ip G (HV)* there exists a unique u G H 1 such that 

(ip,(j))o = / (V« ■ V0 + u<p)ixdm. 
H i Jb 

Formally, the right hand side will be 

/ (V ■ (Vw/i)/!" 1 + u)<fifxdm. 
Jb 

We identify ip as V • (Vti/ij/j -1 + u and the dual pair will be the L 2 inner product. 

Remark 7. With the weight function \i so chosen as f |2.2|) . we observe that if <p G H^, 
then <fi G W 1 ' 1 since 

J (|0| + |V0|)dm < C (J (|0| 2 + |V0| 2 )/ic/m^) ' (J ^ l dn^j ' < oo. 

From the standard trace theorem, the map 

T : flJ(S) — > l}(dB) 
4> !->■ 4>\dn 

o 

is well defined. Thus, the element in H 1 is characterized by the zero trace, and the 
Dirichlet data (I3.3cp makes sense. 

The well-posedness of the w-problem (13 .3p is stated in the following. 

o 

Theorem 8. Suppose that w G L 2 , h G L 2 (if*)* and k G L 2 t with Tt(k) = 0. Then 

o 

the w-problem A3.3\) has a unique weak solution in H such that 

(3.9) IHIw^^CIWIia + IWI'ft,, ) 



with F defined in A2.1$. 



This result when b > 2 and q = was proved in For general case we proceed 
in several steps. 
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An embedding theorem. We define 

r P b / 2 ~\ o < b < 2, 

(3.10) p* = { p- 1 , 6 = 2, 

{ p~ b '\ b > 2. 

We call pf as the conjugate of \i due to the Sobolev inequalities in the following 
lemma. 

o 

Lemma 9. If <ft G H 1 , then 

(3.11) J \(t)\ 2 p*dm <C J (|0| 2 + \V(f)\ 2 )pidm. 
Also, if (j) G for 9 < 1, then for any 5 > 

(3.12) / |0| 2 p~ 1+5 dm < C [ (\<j)\ 2 + \V<j)\ 2 )p e dm. 



Proof. We refer to [25] for a proof of (13. lip when b ^ 2, as well as (13.121) . Here, we 
prove only the case b = 2. 

First for C = maxi< p <2 [pp\~ l we have 

2 /pdm<C [ \<f)\ 2 p J dm+ [ \<P\ 2 /pdm 

B Jl<p<2 J0<p<l 

i r<2 



f f G 

<C \<fi\ 2 p,dm + / — dp 

Jb Jo P 



where we have used the spherical coordinator representation with p = 2 — r 2 and 



(3.13) G 2 {p) = - J im^-'dSr [^y 1 = ±J \4,{rt)\*r»->dSt. 

o 

Note that from <fi e if* one can verify that G(0) = 0. It is known (see [2S])) that 
J g{t)dt^j —dx<C J g 2 (x)x\ \nx\ 2 dx. 

Thus, 

(3.14) f —dp<c[ (G p ) 2 p\ In p\ 2 dp < C [ ^f/idp<C [ {G r ) 2 fidp, 
Jo P Jo Jo r Jo 

where we have used the fact that p\ \np\ 2 < p = pin 2 (e/p). Differentiation of (I3.13P 
in term of r leads to 

2GG r — [ 0V0 • £r N - 2 dSt: + / \<j)(rO\ 2 r N - 3 dS^ 

J\a\=i 2 j m=1 

Squaring both sides and using the Cauchy-Schwartz inequality we obtain 
AG 2 (G r ) 2 < 2 [ <p 2 r N ~ 2 dS^ [ | V0| V^ 2 ^ + {N - 2) ' ( [ <p 2 r N ~ 2 dSt) ' . 
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where we have used the fact r > 1. Hence 

(G r ) 2 < / \V4>(rO\ 2 r N - 2 dS, 4 { ^^G\ 



I5l=i 



which inserted into ( 13.141) ensures that the term —dp is also bounded by C||0||^i. 
The proof is now complete. □ 

Energy estimates. We return now to the bilinear operator B. 

Lemma 10 (Energy estimates). For any t, there exists a constant C which is depen- 
dent on N, b such that 

(1) forw(t,-)eHl 

(3.15) J J \Vw\ 2 pdm < B[w, w; t] 4 C(l + \k\ 2 ) J w 2 pdm; 

(2) foril){t,.)eHl and $ e fy, 

(3.16) »,<MI <C(1+I«I)II^II^II0II^- 
Proof. From (13 .7p it follows 

(3.17) - J Vw ■ Vcpfidm = B[w , (p; t] + J nm ■ V '4>w pdm 4 

where K is given in (13. 6p . 

(1) If < b < 2, then K = 0; hence 

(3.18) ^ |Viu|V*m = B [w, w; t] + J nm ■ Vwwfidm 

< B[w,w;t] + - J |V 'w\ 2 pdm 4 b\n\ 2 J w 2 fid 

and 



Kw(j)dm, 



lm 



\B[i(>,4>;t}\ < \V^j\\V(f)\fidm+Vb\K\ J \i/}\\V<f>\fJidr 

< C(i + |k|)||VH^||V0|U,. 

(2) For b > 2, it suffices to estimate the A'-related term. If b = 2, we have 



K = (N 4 2nm ■ m) In - < IN 4 2&k|)VuM*- 

P 



If 6 > 2, we have 



# = Q-l) p^ /2 (iV4 2 



/tm ■ m) 
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Hence for b > 2 we have 

/ Kw 2 dm < C(l + / w 2 y/ fifi*dm 

This when added upon right side of ( I3.18P using (13.111) with some small e leads to 
( 13. 15H . Using (13. 11 f) again we have 



Kip<f)dm 



<c(i + |«|) J ivil^lv^^^cci + i/cDiHIflill^llfli, 

which when combined with the above estimate for b < 2 gives (I3.16p . □ 
A priori estimate. 

Lemma 11 (A priori estimates). Let w be a weak solution to $3.3\) . Then 
(3-19) sup \\w(t, + < e*W (ikolll, + \\% (hr ) ■ 

with F defined in \2.$, and furthermore 



(3-20) \H\\<e F ^H\\M\ll + \\h\\) 2( - * 



Proof. From the weak solution definition in ( 13. 8 \ we have for any <p G H 1 
(3.21) (d t w 7 <f>)o 1 +B[w,<i>;t] = (h,<f>)o 

By ( I3.16p . (d t w, (f>) o is bounded by 



llH} +C{l + \K\)\\w\\ H ^\\ H}r 

Hence 

(3-22) II^H^). ^ 11^11(11).+ Cf ( 1 + l*DIHIflJ- 

Next we set <fi = w in ( 13 . 2 1 [) and use ( 13 . 1 5[) to have 

M\l* + \ J |Vw|Vm< \\h\\ {ki) J\w\\ Hh + C(1 + |k| 2 )|MI£ji 



1 d „, 
2dt' 



< 2II&II o +-ikn^ +c(i + i«niHi^. 



Hence 

(3-23) + < <7(1 + |«| 2 )||«;||ig + 4],.., 



and therefore by Gronwall's inequality, 
which together with (13T22|) yields ( 13~20|) . □ 



\w \\ 2 L 2 + \\h\\ 2 „ 
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Proof of Theorem We construct a weak solution to ( 13. 3p using the Galerkin ap- 

o 

proximation. Let {0j} be a basis of H 1 and L 2 ^. Then an approximate solution w\ 

i 

in a finite dimensional space is defined as w\ = ^d{(t)&. Here d\(t) is a unique 

i=i 

solution to a system of linear differential equations, 

(d t w h <pj) o + B[w h <pj] t] = (h, <pj) o , 

rf !(0) = ((&, 0i)^(wo, 

where = (0i, • • ■ , 0i) T . Using the same argument as that in the proof of Lemma 
[TIT we obtain estimates for wi such that 

I 2 . . J- ll/wil 2 <r 
'4 (HI 



\ Wl \ \i 2 + \\d tWl \\ z o <e^ |K| M | Kills + ,, 



Extracting a subsequence and passing to the limit give a weak solution w in E. The 
uniqueness follows from the a priori estimate ( I3.20p . □ 

To return to the Fokker-Planck problem (13.11) we will also need the following 

Lemma 12. Let h = —d t q — n~ l L[q\. If q G C]E^ and k G L 2 with Tr(n) = 0, then 

|2 



(3.24) \\ h K* f Sn»^ C (1 + 1*01^)11^. 

L t( H iJ.) Jo 



Proof. For q G G\E\, it is obvious that d t q G L 2 (ifJJ* since E\ C (flj)* C (# r , 

o 

In order to show G L^H 1 )*, we use integration by parts and (13.161) to get 



jj, 1 L[q](pndm 



\B[q^;t]\<C(l + \K\)\\q(t,.)\\ m U\\ Hb G C c °°. 



Since is a dense subset of i/*, for any G i/* with ||0||ij-i = 1, we have 

(3.25) K^Liq]^).] <C(l + |«|)|| g (t,.)||fli. 

Taking the L 2 norm in t leads to the desired estimate. □ 

Theorem [8] and Lemma [T2l lead to the following result for problem (13.11) with general 
Dirichlet boundary condition. 

Theorem 13. Suppose that fou' 1 G L 2 q G C\E X and k G L 2 t with Tr(n) = 0. Then 
for any T > the Fokker-Planck problem h3. 1\) has a unique solution f such that 

(3.26) / = u(w + q) with w G U for < t < T. 
Moreover, for F defined in \2.J$ , 

(3.27) supH^Ct, + < e^CI«D (| |^o| ||, + ^(1 + | | 2 ) | | | 2 d^ . 
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Proof. The estimate ( I3.27P follows from (I3.19P and the estimate in Lemma [T2| with 
Fe F replaced by e F . 

We now prove uniqueness of / in terms of q\oB- Let fi(i = 1,2) be two solutions 
with qi such that q\\dB = <?2 \bb an d same initial data / . Set w = (/ 2 — fijv" 1 , then 
u> solves w-problem (13. 3p with «; = /i = 0. Hence w = 0, leading to /i = / 2 . □ 

o 

Remark 14. yls mentioned in Section^ that H 1 = if b > 6,i.e., the trace of 
q G i?* vanishes if b > 6. TTras, i/ie boundary condition f)3.1cl) zs nothing but a zero 
Dirichlet boundary condition. In Section® we show the well-posedness with a nonzero 
Dirichlet boundary condition for b > 6 using yet a different transformation. 

3.2. Probability density function. So far we have discussed well-posedness of the 
initial-boundary value problem fl3.ll) for b > and any given q. We now turn to the 
question of which q corresponds to the probability density, i.e., non-negative solution 
with constant mass for all time. 

Proposition 15. Let f(t,m) be the solution to problem (13. ip obtained in Theorem 
[TR If fo > and q(t,m)\dB > almost everywhere, then f remains nonnegative for 
t > 0. 

Proof. We adapt an idea from [7]. Let / be the positive and negative parts of the 
solution / such that f = f + — f~ ■ Obviously, w ± := f ± i'~ 1 G H'h and q\ss > 0. This 
implies that the trace of w~ at the boundary vanishes, so 

w~ G Hi. 

From the equation 

dtw/i + L[w] = 0, 
which is transformed from (I3.1al) it follows that 

(d t w, w~) o + B[w, w~; t] = 0. 
"it 

Since (dtw + ,w~) o ^ and J L[w + ]w~dm vanish, hence 

"it 

2d~i\J \ w ~\ 2 f Ji ^ m ^j ~^ ^[ w ~ ^ w ' A = ®- 
The coercivity of B, (13 . 1 5j) . gives 

Hence 

sup|K(t,-)||i, <\\wo\\he F(]K]) 
t M M 

for T > 0. Since Wq = 0, \\w-(t, -)\\ 2 L2 = for all < t < T. □ 

Proposition 16. Let f be a solution to the Fokker-Planck problem (13. ip obtained in 
Theorem\]J^ Suppose b>2 and q{t,m)\dB > 0. If q\dB = for all t G [0, T], then 



f(t,-)dm= fodm, te[0,T], 



and vice versa. 
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Proof. It suffices to prove the claim for smooth enough / since the general case can 
be treated by an approximation as in We rewrite ( 13. lap as 



dtf = -V ■ (Kmf) + V ■ (p b/2 Vjj^J 



First, we take a test function 4> £ (m) = <f> e (\m\) G C%°(M. N ) converging to \b as e — > 
such that 

I 0, |m| > Vb-e/2 ' 1 re| ~ e 

and for any smooth g 

p Vb-e/2 

(3.28) / g(r)(j)' £ (r)dr -»■ as e 0, 

where 0' (r) = V</> £ • , — r. 

|m| 

One can construct such a £ by mollifiers, for example 

4> E {m) = I Ve/^( m ~ m')dm! 

where 



B Vb-3e/4 



rj e (m) = -^rj(m/e), rj{m) 



f Ce HH a , 


\m\ 




\m\ 



and C is the normalizing constant. 

Since V0 e is supported in B e := B^_ £ , 2 \ -^%/&-e> nence 

(3.29) j t j f<p e dm = j /«m ■ V<Mm - ^ P V2 V f • V0 B dm. 
By ?y = /f -1 , the right hand side reduces to 

(3.30) / (wKm - Vw) ■ V(j) £ vdm - I wp b/2 V(f) £ ■ V(up~ b/2 )dm. 
Jb e Jb s 

The first term converges to 0. Indeed, 

j \wKm — Vw\ 2 fidrn^j (^J |V0 e | 2 — drn^j 



{wnm — Vw) ■ Vcf) £ vdm 



< 



Since v 2 / p, = p b l 2 for b > 2, by mean value theorem there exists r G (y/b — e, y/b — e/2) 
such that 

/ \V^ £ \ 2 -dm = £ - [ \V^\ 2 p b ' 2 dS<CE h l 2 -\ 

JB e P 2 J dBr 

which is uniformly bounded for b > 2. Using w G H^, we obtain f BS \wnm — 
Vw\ 2 fidm y as s y . Hence the first term in (13.301) converges to 0. 
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On the other hand, for Co 



-2, 6 = 2 
2-6, 6 > 2 



L 



wp b/2 V(j) £ ■ \7{up~ b,2 )dm = C wW<p £ ■ mdm 

J B? 



pVb-e/2 r 

Co / wr(f)' £ (r)dSdr 

J -Jb-e JdB r 
rVb-e/2 / r \ 

Co I [r wdS j <p' £ {r)dr. 



Vb-e \ JdB 



Due to (I3.28P this converges to 

-C Vb [ wdS = -CoVb [ qdS. 

JdB JdB 



Since C ^ 0, this shows that — / fdm = if and only if / qdS = 0, or q 

' B JdB 



dt 



<1B 



0. □ 

Remark 17. In Proposition fig), the assumption b > 2 is sharp. In the case b < 2, 
we need to consider nontrivial q ^ since the equilibrium profile f eq = p h l 2 satisfies 

q\ dB = p^ 1 \ dB = l. 

This requirement is also consistent with [35], in which it was shown that when 6 < 2, 
f l>1 \dB — q\dB is necessarily prescribed and each solution depends on the choice ofq. 
It would be interesting to figure out a particular q for which the corresponding solution 
when b < 2 is a probability density. 

4. The Fokker-Planck equation 

In this section, we show the well-posedness of the FPE (12.91) including x variable. 
The result is stated as follows. 

Theorem 18. Suppose that for 6 > and any integer s > N/2 + 1, V ■ v = and 

(4.1) v e C t H s x n L*H?\ fo^ 1 e H S X L% q e C l t Hl +1 Hl 0<t<T 
for any T > 0. Then ( \2.9\\ has a unique solution f = v[w + q) satisfying 

(4.2) sup ML + \ f \V m w\l,sdt < e F(Hs+l) (Nl^ + IMI? lW ) > 

t * Jo 

where F was defined in {2.1$ . 

The proof of Theorem [18] consists of two parts: first we show the existence of the 
solution / to problem (12. 9 p by using the flow map, followed by proving regularity in 
x inductively such that w G CtH^L 2 ^ n LfH^H^ with v, fo and q given in (14. ip . In the 
second step, we derive estimate (14. 2p directly from (12. 9p to control / in terms of the 
given data. The uniqueness can be obtained from the estimation (14.21) as performed 
in the proof of Theorem [TS] 

First, we state a technical lemma. 
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Lemma 19. Suppose that if; £ and <fi £ H 1 . Then for the trace map T : 
W l >\B) -> L\dB) 

(4.3) T(^0/i) = 0. 

o 

Proof. Since is a dense subset of H 1 , it suffices to show that for a fixed ^ £ if* 
and any £ C£° 



(4.4) ||^|| iyl , 1 <C'||0|| Jf i. 

Then, the standard trace theorem asserts that T(ip<frfi) is well-defined in L 1 (dB) and 
it vanishes, also T is a continuous map with respect to <p, we can thus conclude (14. 3[) 

o 

for any (f> £ H 1 by passing to the limit of sequence (p n £ C£° such that <f) n —> (j). 

(I4.4p is indeed the case. It is obvious that V m ?/;^ and ^V m 0/x are integrable. 
For 6^2, |V m /i| < CVw? and d3~TT]) yield 



For 6 = 2, 



|^0V m Atdrri| < C||V|U?J|0||jyi- 



|V m /i| < C(ln 2 - + In -) < C(ln 2 - + y/jljF). 
P P P 



Using (I3.12p and ip £ if* we obtain tp £ L?_ 1+s for any 5 > 0. Hence 



^0 In — dm 



~°yj I ^' 2p ~ 1+5dm \l / 0V _5ln4 (^)^ 

It follows that for any 6 > 

1^1 + |V m (V>#)|dm < C||V||^||0||ifi 
as we desired. □ 

The main ingredient for the proof of Theorem [18] is to use the calculus inequalities 
in the Sobolev spaces, see Appendix 3.5 of [3D]: for any positive integer r > and 

(4.5) 2^ Wd^iuv) — ud J v \\ L 2 < C (\\Vu\\L°°\\v\\Hr-i + \\u\\ht\\v\\l°°) , 

h\<r 

(4.6) ||wv||//r- < C(j jii| |z,«>] ]u] l/r^ + I M l-ff r l M 

Note that (14.51) remains valid when d 1 on the left hand is replaced by the correspond- 
ing difference operator. 

Proof of Theorem Q3 

Stepl (well-posedness) Let a particle path be defined by 

d t x(t,y) = v(t,x(t,y)), x(0,y) = y, 
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along which the distribution function f(t,y,m) := f(t,x(t,y),m) solves 
(4.7a) dJ+Ctf] = 0, 

(4.7b) f{0,y,m) = fo{y,m), 

(4.7c) f(t,y,m)v- l \ aB = q(t,y,m)\ 9B . 



Here £ is defined in (13. 2p with n replaced by n(t, y) = Vv(t, x(t, y)), and q(t, y, m) : = 
q(t,x(t,y),m). 

In order to show existence of the solution to (12. 9p under the conditions v G CtH s x H 
L 2 H^ +1 and V • v — 0, it suffices to show that (14. 7[) has a solution / = u(w + q) such 
that 

w := w(t, x(t, y),m) G C t HyL 2 ^ fl L 2 H y H x , 

assuming that 

(4.8) k G L 2 H y) w eH s y Ll qeClHffi. 

These follow from (14. ip since < C\v(t)\ s+ i for t > 0, ioo(x, m) = /o^ -1 — — 

0) = Wo(y,m), and ||?||i,i, s < C||g|| i 5 i ;S _|_i, for which we have used d t q = d t q + v ■ Vq. 
Using Theorem [13] for each y, there exists a unique solution / such that 

f = v{w + q) 

with w satisfying ( 13. 2 7ft . i.e., 

sup | \w(t, y, •) + l IK-, •) ||i f m < e 15 '^!) (| My, 

(4.9) +^ T (l + |K(-, 2/ )| 2 )||g(t, Z /,-)|| 2 ^. 

Integration of (14. 9 j) with respect to y, upon exchanging the order of integration in y 
and m, and using the Sobolev inequality, sup y |k| < C|/c| s _i, gives 

(4.10) sup \w\l fl + \ ! \w\\ fi dt < ^(N-O (| Wo |2 + \\q\\l lQ ) . 

Hence w G C t L 2 L 2 fl L 2 L 2 H^. On the other hand, the right hand side of (14.91) is 
uniformly bounded in y, taking sup y of ( 14.91) gives 

(4.11) sup \\w(t,y, < e^-^kolL-i + 1 1^1 l?.x..-i) - 
t,y " 

We now use an induction argument to prove that w G C t H r y L 2 S\l^H T y H x for < r < s, 
and 

(4.12) sup H; j +i [ |tS|? ir dt<e F ^(|«;o|g 1# + ||?||? il , a ). 

t 2 J 

The case r = has been proved as shown in (14.101) . Suppose (14.121) holds for r = k, 
we only need to show (14.121) for r = k + 1 < s. 

To prove regularity of / in the y variable, we use difference quotients. Define the 
difference operator in the y variable as 
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Apply <5 7 to (H2D with \-f\ < s, then 

(4.13a) dtFf + C[iPf\ = V m -J, 

(4.13b) 5 7 /(0,2/,m) = 5 7 /o(2/,m), 

(4.13c) 5 7 /(f,?/,m)z/ -1 | 9 B = S J q(t,y,m)\ dB , 

where 

(4.14) J = m5 7 / - Ffcmf). 

This when transformed into the w-problem of form (13. 3p involves the following non- 
homogeneous term 

(4.15) h = -d^q - H^L^q] + V m • Jv' 1 . 

Using Theorem [T3] again for each y, 5"' f is the unique solution to (I4.13P as long as 
h E L 2 t (H^)*. Moreover, 

where 5^w, using (I3.19p . satisfies 

sup \\Fw(t,y, -Mi + h\SM-:V, -)\\% H i < e Filk[ -' y)l) (ll^oWh + IM** 
Integration in y gives 



sup |5 7 u;| 



T 



(4.16) 



\^w\l dt < e F (™Py\~<;y)\) ( |^ |o,o + \\h 



We now turn to bound the last term in the above inequality. For any <p G H 1 and J 
defined in (I4.14p . Lemma [T91 allows the use of integration by parts. Hence, 

J V m • Jv~ l 4>iidm < | Jis~ 1 \\isV m — \\<f>\dm + J | Jz/ _1 | | V m 0|/i<im^ 

^ciiJ^l^di^ii^ + nv^ii^) 
^ciiji/^ii^ihIbi. 

Here we have used |^V m — | < C\f\f\x and the embedding theorem (13. lip . This 
together with Lemma [T21 and (I4.15P yields 

(4 ' 17) 

\\hf „ <C f (l + supl^y)] 2 ) f\\Pq(t,y,-)\\ 2 dydt + C [ | Ji/" 1 ]g dt. 

L t L yi H ]i)* Jo y J Jo 

For |7| < s, the first term on the right side is bounded by 



(4.18) 



n\*\s-iW<i\\Uo <^u)ikiii 
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To obtain (I4.12p for r = k + 1 < s, it remains to estimate the last term in ( I4.17P with 
I7I — k + 1. In fact, 



\Jy- l \U = \{5\Rmf)-km5y)u-X fl 

< C(sup |V y «| 2 |/i/ -1 |g jfc + sup H/^l!^ 

y y 

< C\k\ 2 s (\w\ 2 k + sup \\w\\ 2 L 2 + ||g||i,i, a ), 

y 

where we have used (14.51) with <9 7 replaced by <5 7 . 
Using (14.121) for r = k and (I4.1ip we have 



ri^- 1 llo*<e Fa * W (Kll. + ll?llli 1 .). 
Jo 



This and 04.181) when inserted into (14.171) gives a bound for II h\\ 2 . That bound 

L^LliHl)* 

combined with (14.161) yields 

sup l^tSlJo + \ f \Pw\\ fi dt < e F ^ (\w \l s + \\q\\l 1>s ) < 00, | 7 | = * + 1. 
t ^ Jo 

Sending rj — > we obtain (I4.12p with r = k + 1. Hence, (14 . 1 2 j) holds for any r < s, 
and thus the solution / to ( 12. 9 p exists, and 

1 f T 

SU V\ W \l,s + o / IHl,*^ < °°' 

t * Jo 

One may obtain an upper bound from (I4.12p with r = s using the inverse map of 

x = x(t,y). Nevertheless, the next step gives the claimed bound in (14. 2p . 

Step2 (a priori estimate) For a priori estimate, we consider the w-problem (12. lip 

(4.19) fi(d t + v ■ V)w + L[w] = -n(d t + v- V)q-L[q]. 
Recall that 

L \ w \ = ~\^ m ■ (V m w/i) + V m • (Kmwfi) - Kw. 
Take 7 derivative in x-variable. Then, the left and right hand side of (14.191) will be 

(4.20) / = /i(d t + v ■ V)d 7 u> - ^V m ■ (V^ai/i) 

(4.21) + ^((v ■ V)w) - ■ V)d^w] 

(4.22) + V m • (<9 7 («;mu;/i)) 

(4.23) - d^(Kw), 

(4.24) // = -nd t 8T'q+'^7 m .{V m 8rqn) 

(4.25) - fid"<((v ■ V)q) 

(4.26) - V m • (d^nmqfi)) 

(4.27) + d\Kq). 



21 

We now estimate term by term of 



(4.28) J J Id^wdmdx = j j ndlwdmdx - 

M<« \i\<s 

Since v is divergence free, the first two terms on the left hand side will be 

1 d . ,o 1 ,_ l2 
2^Mo, s +2|V^| 0l , 

Indeed, Cauchy inequality shows that the term related to (14.21 j) is bounded by 
e\w\l s + C £ J2 [ f \d 1 {{v-V)w)-{v-V)d 1 w\ 2 ^dmdx 

Now, we exchange the order of integration in x and m, and apply ( 14. 5 p to obtain 
e\w\o,8 + C e J (l|Vw||! ? ||Vw(-, + \\v\\ 2 H s\\Vw(-,m)\\ 2 L J) fidm 

^ I 1 2 i I 1 2 1 1 2 

where the Sobolev inequality, |w|o < C|w| s _i for any u G H* -1 , is invoked in the last 
inequality. Similarly, the term with (I4.22p will be estimated as follows due to (14. 6P ; 

\l\<8 



< s\V m w\l iS + C £ J (\K\\^\w(-,m)\l + \K\l\w(-,m)\\^) iidm 

< e\V m w\ 2 0tS + C e \v\ 2 s+1 \w\l )S 

Recall that 

{0, 0<b<2, 

(N + 2ktti ■ m) In-, b = 2, 
P 

(N + 2Km-m)(b/2- l)p l ~ b/2 , b>2. 

Thus, we can express K as 

(4.29) K — ciV ' im* + c 2 Km • m^fifi* 

for some positive constat q depending on N and b. We now estimate the last term 
on the left hand side, by using 

d 1 (Kw)d 1 w = Cild^w^y/JIJF + C2<9 7 (mti • mw)d' y wy/njj*. 

The Cauchy inequality and the embedding theorem A3. 11 j) give 

c = * / w«, m )i»™ 

< £ / Wt ,,»)iyjm+C,/w,-,m)iy™ 
< e|V m iy|n s + C e |w|o 
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Similarly, 



The last term, using (14. 6 p and the Sobolev inequality for k 

by 



Id' 1 (him ■ mw)\ 2 [idmdx . 



Vi>, is then bounded 



C e\v\ 2 s+l \w\l s . 



Hence, 



/ / d 1 {Kw)d 1 wdmdx 

\l\<* 



<£|v OT ^lL + c e (kl' + i + i)hlL 



Now we turn to the right hand side, related to (I4.24I) - (I4.27I) . The estimation is 
similar to that for the left hand side. Except that here we have to assume higher 

regularity of q in x than that of w since / v ■ / V 'd 1 qd 1 'w [idmdx does not vanish as 



v ■ / V ' d 1 wd 1 w [idmdx . Indeed, the first two terms, related to (I4.24p are bounded 



by 

+ c £ \ 

and the other terms are estimated as follows; 



h\<s 



d 1 {v ■ W q)d 1 w [idmdx 

/ / d 1 {i\mq)V md 1 w [idmdx 
7 |<, J J 

y y d 1 {Kq)d 1 wdmdx 



^ Yi i /^f I |2| |2 

< eHo, 8 + c 'eM.Mo,H-i 



< e|V m w|o. + C e |v|2 +1 |g 



0,s> 



< e|V m iu|o. + C e |5|L + C e |i;|2 + i|g| 



We combine all estimates for sufficiently small e to obtain 

(4.30) d t \w\l,s + l\V m w\l s < C(\v\ 2 s+1 + 1) (ML + + 



We deduce that 

|2 



\W\ 



l \v m w\l s dt < ^ (holt + • 

Replacing Fe F by e F leads to (14. 2p . 

5. Coupled system 



□ 



In this section, we prove Theorem [T] by the fixed point argument as described in 
Section [2J 

We begin with a key lemma, which will be used to estimate the stress r. 
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Lemma 20. Suppose that <p G H 1 . For any e > there exists C £ such that 



(5.1) 



4>up l dm 



<a 



fidm + E / |V m4>\ 2 ^dm. 



Proof. For b > 2, the Cauchy-Schwartz inequality yields 



tm 



< J \(j)\ 2 ^dm J /i 1 dm. 



For any e > 0, taking C. = /„-*» < oc, we obtain (ED for ft > 2. 



For 6 < 2, we define for fixed M, 



G = {0 G Hi : 



~ l dm 



HI < M}. 



It suffices to prove 



I := inf 

4>eG 



'jjidm > 0. 



Let {4> n } C G be a sequence such that 



lim 



= inf 

0eG 



2 fidm. 



Since {</>„,} is bounded in if*, by embedding theorem (13.111) . there exists a subsequence 
{<pn k } such that 



Furthermore, since 




'"ft 



G Lt for b < 2 

A 1 



in # 



in LJ.. 



y 0Vp 1 dm = J ( l ) *\p^> 



—[i dm 



= lim 

This shows that 0* G G. On the other hand, 




/*l 2 



jidm < lim / |0 n J 2 /x<im = /. 

n fe ->oo J 

If / = 0, then 0* = which is a contradiction to 0* G G. 



□ 



The zero trace of is essential for the estimate (15. ip . For the general case, i.e., for 
G H^, one can only have a weaker estimate. 

Lemma 21. If <fi G , £/ien t/iere exists C such that 



(5.2) 



< C 



HI' 
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Proof. For b > 2, we have 



J <\>vp 1 dm < C J \(p\ 2 pdm, C := J p 1 dm < oo. 



For b < 2, 



' x dm 



<c* 



l2 p- 1+s dm, C s :-- 



u 2 p- l - 5 dm 1 . 



We choose S > small enough so that Cs is bounded. On the other hand, by (13.121) 
in Lemma [9] we have 



J \<f)\ 2 p- 1+5 dm < C J{\<f)\ 2 + \V m <P\ 2 )p b/2 dm = C 

2 + \V m 4>\ 2 )pdm < C 



i)\ 2 p~ l+& dm < C 

This completes the proof. 
We now turn to the map 

and 



2 + \V m (f)\ 2 )pdm, b<2 
2 + |V m 0| 2 )pdm, 6 = 2. 



□ 



T : M -> M 

(m, -Cct) i—t- (u, wl 



M={(v,iw): sup < Ai, sup ML + « / I V m iu|jj s dt < A 2 

I 0<t<T 0<t<T <t Jo 

We first prove that, given v G H s x , j§v~ x G H^L 2 ^ and q G C t 1 if| +1 i^, the map J 7 is 
well defined, i.e., .F(M) C M for some A 1 ,A 2 ,T. 

Let (w, cut) G M. It is now well known that ( 12.81) has a unique solution v such that 



T 



(5.3) sup|v|*+ / MLi<J* < \v \ s + C / |u|>|*dt+ / \r\ 2 dt, s > N/2 + 1 



By Gronwall's inequality and sup \u\ s < Ai, we have 

0<t<T 



(5.4) 



sup \v 

t 



v\ 2 s+l dt< \v Q \l + 



\r\At I e 



We proceed to estimate the stress term 



f \r\ 2 s dt= [ [ \d J T\ 2 dxdt 

Jo Jo |7| < s 7 



where using Lemma 



m ® md J (w + q)vp l dm 



<CeJ \& 1 w\ 2 pdm+ E - J \d^V m w\ 2 pdm + 2b 4 J d^qvp^dm 
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Using (15. 2 p the last term is uniformly bounded by 

c\\dMt,x,-)\\h,< ciMlli ,s+l- 

Hence for (w, ro) 6 M we obtain 



(5.5) 



f \r\ 2 s dt < C £ TA 2 + eA 2 + CT|g|? ilf-+1 < CT(A 2 + ||g||? iW ) + eA 2 , 
Jo 



where we have used the assumption q G C^H^ 1 !!*. 
We choose Ai as 

(5.6) A 1 = 2\v Q \ 2 8 e, 

A 2 as 

(5-7) A 2 = (\w \l s + \\q\\l liS+ Ae c ^ 

for T < l/(Cy/jQ. 

Hence, if T and £ are chosen small enough so that 

CT(A 2 + \\ q \\l hs+1 ) + eA 2 < j-Ax, 

we get 



(5-8) e c ^ T {\v \l + CT(A 2 + |g|? )1)S+1 ) + eA 2 ) < e(\v \ 2 s + -Ax) < ^. 



^ T {\v \ 2 + CT(A 2 + l-l 2 ^^U,J2^1 
This together with (15. 4ft . (15. 5 p gives 
(5-9) supH^+ / \v\ 2 s+1 dt < A 1 . 



Estimate in Theorem [HI (JEZD and (ESD yield 



(5.10) 



1 f T 

su PklL + o / |V m w| dt<^ 
/ 2 J 



So the map J 7 is well defined in M. 

Next, we show that J 7 is a contraction mapping for small enough T using a weak 
norm on M, i.e. 

1 f T 

(5.11) ||(u,w)||^:=sup|u|g + supHj| +- / \V m w\ 2 0fi dt. 

t t £ Jo 

Suppose that Vi(i = 1,2) are solutions of the NSE (12. 8 p with Ui(i = 1,2) and Tj(i = 
1,2) computed from Wi(i = 1,2) respectively. Then we obtain 

(5.12) d t v + {u 2 ■ V)v + {u - V)«i + Vp = V ■ r + Av, v(0, ■) = 0, 

where v = v 2 —v±, u = u 2 —u\,p = p 2 —p\, r = t 2 — t\ and w = w 2 —w\. Multiplication 
by v to (15. 121) and integration with respect to x yield 

- — Mo + J (u ■ Vv i)vdx = — J rVvdx — J \Vv\ 2 dx. 
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(5.13) 
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^Mo + |V^ < \u\l + |r|g + sup|Vr;i| 2 Hg 



Let fi be the solutions to (12.91) associated with Vi{i = 1,2). Then 

w = (f 2 - =: w 2 - Wi 

solves 



(5.14a) d t wfi + v 2 ■ Vu>yU + L 2 [w] = —v ■ Vu>i/i — V m • (Vvmwiis) 

(5.14b) w(0,x,m) = 0, 

(5.14c) w(t, x, m)\gB — 0, 

where L 2 [w] = defined in (13. 5p with k = Vv 2 . Note that Wi\gB = <?|aB, i-e. 

o 

Wi(t, x, ■) G H*, so w(t x •) G 
We deduce from (I5.14aj) that 

- — \w\q o + -| V m u>|g o < J J \\7v 2 m ■ X7 m wvw\dmdx + J J \Kw 2 \dmdx 

+ y ' |i> • y Vwiw/j,dm\dx + ^ ^ V m • (Vvmwiis)—wdm 

Similar to that led to (14. 2p . first two terms on the right hand side are bounded by 

C £ (\v 2 \ 2 s + l)\w\l + e\V m w\l , 



dx. 



and the third term 



< C|«|>,|^ + | m |? |0 . 

The last term, using integration by parts with vanished boundary term due to Lemma 
PT9l is bounded by 



J J V m • (yvmwiv)—wdm dx = J J S/vmwiV ■ V m {j~ w ^j dm 



dx 



< C e \Vv\l\ Wl \l s + e\V m w\l . 



Putting all together we have 
d 



-rJ W \o,0+n\V m w\Z Q < C(\v 2 \l + l)\w\l + C\ Wl \iMt+\Vv\l) 



dt 



< C(A 1 + l)\w\ 2 00 + CA 2 (\v\ 2 +\Vv\l). 



Substitution of the estimates of \Vv\q and j^\v \% in (15.131) gives 

(5.15) -(\v\l + Ho, ) + l\V m w\l < D(\v\l + \w\l ) + D\u\l + D\rg 

where D is a large constant depending on C, A\,A 2 , for example we may choose 

D = C(A l + l)(A 2 + l). 
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The Gronwall inequality gives 

sup(|i;|jj + Moo) + I [ \V m w\ 2 dt < De DT [ \u\ 2 + \r\ 2 dt 
t * Jo Jo 

for any < T* < T. Due to the similar estimate for r as (15. 5p . the right hand side is 
bounded by 

De DT I T* sup |u|o + C £ T* sup \w\q + e / |V m 07|o di I . 
V t t ' Jo 'J 

We choose e = - J~ — , T* = ^min<T, — \ 1 and redefine T = T* to 

ADe DT 2 [ (C £ + l)De DT ) 

obtain 

(5.16) ||(va,u7 2 ) - («i,u;i)||m = II(^ w )IIm < TjllK,^) - («i,^i)||m. 

This shows that J 7 has a fixed point (v, w) in M, which is a solution to the coupled 
problem (12.11) . Since jF(y,w) = (v,w), (15. 3p and Theorem [131 imply that (v,w) G X M . 

The uniqueness follows from the same computation of estimates for the contraction 
mapping. Let {v ir /j^ _1 )(z = 1,2) be solutions of the coupled problem (12.11) . Then 
v — V2 — v i solves (I5.12p with Ui = u = v, and r = r 2 — T\ computed from fa. 
w = (/ 2 — fi)v~ l also solves (15.141) with = j\v~ x . Similar to (I5.15p . we obtain 

"^(Mo + IHo.o) + \\VmW\l,0 < D (\ V \l + IHo,0 + Mo)- 

It follows from the estimate for r and Gronwall inequality that (v, w) = (0, 0), which 
gives the uniqueness of problem (12. ip . 

6. A FURTHER LOOK AT b > 6 

In this section, we sketch proofs of Theorem[2]and TheoremH]for the case of p = po- 
Consider (12. 9p when x is not involved, i.e., (13. ip . The corresponding w-problem for 
w = fis" 1 — q with p — no solves (13. 3 j) with the operator L replaced by 

(6.1) LqH = -~V ■ (ywfi ) + - - 6 J m ■ Vwp 9 ' 1 + V • (nmwpo) ~ K w, 
where 

(6.2) Ko = [N(b/2 - 1) + 2nm ■ m(l - 6)} p 6 ' 1 . 

Define the conjugate of po a s (I3.10p . /ig = p ~ 2 , then .ffo can be rewritten as 



(6.3) Kq = [N(b/2 - 1) + 2«m • m(l - 0)] ^0/4 

To ensure well-posedness of (13. 3p . we need to check the coercivity of B [w, w; t], 
which is defined as 

- J \Vw\ 2 podm = Bo[w,w]t] — ^2 — -b — J m -Vwwp e ~ l dm 

V • (nrnw p )wdm + J K w 2 dm. 
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From the proof of Lemma [TOj the last two terms are bounded by 

C £ / w 2 p dm + e / \Vw\ 2 [i dm, 



where the embedding theorem (13.111) has been used. For small enough e, this estimate 
yields 



as long as 



- J \Vw\ 2 p^dm < Bo[w,w;t] + C J w 2 podm, 
J \ 2- h)-e\m- Vwwp^dm > 0, 



for w G H,. n . This is indeed the CclSG, clS shown below. 



' AH) 



Lemma 22. Let w G H},. Then 

H'O 



(6.4) J (2 - h) - 9)m ■ Vwwp e ^dm > 0. 

Proof. From — 1 < 9 < 1 and b > 6, we see that (2 — 6/2 — 6) < 0. It suffices to show 

J m ■ \7wwp e ~ 1 dm = — J m ■ Vw 2 p e ~ 1 dm < 0. 

Integration by parts gives 

/ m ■ Vw 2 p e - 1 dm = - w 2 (Np e - 1 + 2(1 - 9)\m\ 2 p e ~ 2 )dm + / w 2 p 6 ~ 1 m ■ ^-dS 
J J JdB H 

< Vb [ iv 2 p e - 1 dS = 0. 

JdB 



<dB 

Here we use the fact that w 2 p e ~ 1 G W 1 ' 1 and w 2 p e ~ 1 \dB = 0. To see this, for any 

o 

w G H}.„, we estimate 



■ MO' 



J w 2 p 6 ' 1 + IV^V" 1 )!^ < J w 2 p e ~ 1 + 2\wVw\p 6 - 1 + 2(1 - 9)\mw 2 \p e - 2 dm 



< 



C J w 2 ^pqPq + \w\\Vw\y/popQ + w 2 p* Q dm 



< C\\w\\ 2 



H 1 ' 
"mo 



due to the embedding theorem (13.111) . Thus w 2 p e 1 \gs G L 1 {dB) from the trace 

o 

theorem and it is zero from the fact that is a dense subset of H 1 . Thus (16.41) 
follows. □ 

We now turn to the FPE problem including x-variable. The first step in the proof 
of Theorem [T8l remains valid for p = po- To check the second part of the proof, we 
need only look at two extra terms beyond those in 04.28)) . 

2 - h - 9 J I m-V m d J wp e - 1 d y wdm, - (2 - ho - 9] I m-V^qp^d^wdm. 
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The first term is non-positive from Lemma [22j and the second term is bounded by 



C 



J J m-Wmd^qp 8 l ff 1 wdm <C £ J |V md 1 q\ 2 podm + e J ^w^p^dm. 

These ensure the same estimate (I4.30p and thus (14.21) . 

For the well-posedness for the coupled problem, we utilize 9 < 1 and Lemma 
For example, for the proof of Lemma [201 with po 

l dm\ 2 = | / (pdm\ 2 < I (p 2 p dm \ p^dm. 



Since 6 < 1 we have J pQ l dm < oo, hence (15. ip . Verification of other terms is omitted. 

The remaining is to show Theorem HI the solution / is a probability distribution if 
and only if q\gs — for p — p , Positivity of / follows as in Proposition [151 For the 
conservation of mass, as in Proposition ITS} we only have to check (13.30 p . 

/ [wkvii — Vw ■ W(f) e v)dm — I wp b / 2 'V4 > s ■ V(i>p~ b / 2 )dm. 

J B E J B E 

Since u 2 /fi = p 2 ~ e and 2 - 9 > 1 

£ - [ \vj> E \ 2 P 2 - 9 ds 

1 J a B r 

converges to as e — > 0. Thus the first term converges to as well. On the other 



hand, the same argument shows that the second term converges to C / qdS for 

JOB 

some nonzero constant C. Hence, we conclude Theorem [4] under the assumption of 
Theorem [2j 

7. Conclusion 

In this paper, we have analyzed the FENE Dumbbell model which is of bead- 
spring type Navier-Stokes-Fokker-Planck models for dilute polymeric fluids, with our 
focus on developing a local well-posedness theory subject to a class of Dirichlet-type 
boundary conditions 

fv~ l = q on dB 

for the polymer distribution /, where v depends on b > through the distance 
function, and q is a given smooth function measuring the relative ratio of f jv near 
boundary. We have thus identified a sharp Dirichlet-type boundary requirement for 
each b > 0, while the sharpness of the boundary requirement is a consequence of the 
existence result for each specification of the boundary behavior. It has been shown 
that the probability density governed by the Fokker-Planck equation approaches zero 
near boundary, necessarily faster than the distance function d for b > 2, faster than 
d\lnd\ for 6 = 2, and as fast as d h l 2 for < b < 2. Moreover, the sharp boundary 
requirement for b > 2 is also sufficient for the distribution to remain a probability 
density. 
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